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ABSTRACT 

We describe an emission model for gamma ray pulsars based on curvature 
radiation-reaction limited charges in the outer magnetosphere. We show how pair 
production on thermal surface flux can limit the acceleration zones; estimates for the 
efficiency of GeV photon production and the 7 -ray beaming fraction are derived, 
including their dependence on pulsar parameters. In general t / 7 increases with pulsar 
age, but is decreased for low magnetic fields and for small magnetic inclinations. 

We argue that this produces GeV pulse profiles, curvature spectra and detection 
statistics consistent with the observations. We also describe the optical through 
X-ray pulsar synchrotron spectrum and the spectral variations with pulsar phase. 

A test computation for Vela-like parameters reproduces phase-resolved GeV spectra 
consistent with those observed by EGRET . Finally we comment on very high energy 
pulsed emission and particle production and note extensions needed to allow a more 
complete pulsar model. 

Subject headings: pulsars — gamma rays — radiation processes 


1. Introduction 


Isolated spin-powered pulsars are very prominent in the high energy sky. In the GeV 7 -ray 
range, in particular, the EGRET experiment on CGRO has now detected pulsed signals from at 
least seven objects. In addition ~30 point sources are found at low latitudes (b< 10°) along the 
galactic plane (Catalog 2, Thompson, ef al. 1996). The identification of one of these objects 
(Geminga=PSR 0630+18) with a pulsar not visible in the radio has made the connection with 
young isolated pulsars particularly appealing. 

In the outer magnetosphere picture of Romani and Yadigaroglu (1995; hereafter RY), we 
developed a model for the beaming of high energy 7 -ray emission that reproduces the observed 
properties of individual 7 -ray pulsars, including high energy pulse profiles, the relation to the 
radio pulses, and the non-detection of pulsars such as PSR B0656+14. In Yadigaroglu and 
Romani (1995; hereafter YR), it was shown that the statistics of 7 -ray pulsars and Geminga-like 
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objects were in good agreement with the predictions of our 7-ray pulsar model for a pulsar 
properties entirely consistent with those inferred from radio surveys. These results extended earlier 
population studies based on extrapolation from the observed pulsars (cf. Halpern and Ruderman 
1993, Helfand 1994). However the luminosity law used in the population estimates of YR was 
based on a standard phenomenological fit of efficiency vs. spindown age; a more physically based 
description of the gamma ray efficiency is needed to follow dependence on magnetic field and other 
pulsar parameters. 

A second class of important pulsars models are those placing the emission zone at the polar 
cap, with emission altitudes of ~ 10 4 cm to several neutron star radii (Harding and Daugherty 
1993, Dermer and Sturner 1994). While these polar cap models have difficulty in reconciling the 
observed 7-ray pulses with the known radio pulsar geometry, and while substantial modification of 
the standard polar cap picture is needed to derive reasonable pulse profiles and pulsar numbers, 
these computations have had good success in obtaining high energy spectra similar to those 
observed for the Vela pulsar and other young objects. Moreover both the resonant Compton 
scattering based model of Sturner, Dermer and Michel (1995) and, especially, the curvature 
radiation based model of Daugherty and Harding (1996) have shown that polar cap radiation 
schemes can produced spectral index variations through the observed pulse similar to those 
measured for Vela by EGRET (Kanbach et al. 1994). Similar results have been difficult to obtain 
for radiation from the outer magnetosphere. Chiang and Romani (1994) were able to show that 
because of the variation of the magnetic field through the radiation zone, spectral index changes 
were expected during the pulse in a radiation model based on the outer gap picture of Cheng, Ho 
and Ruderman (1986; hereafter CHR). However, they also found that this scheme failed to produce 
self-consistently the soft photon fields needed to support the gap closure, even for the Crab pulsar, 
so that further study of radiation processes in the outer magnetosphere were required. 

In this paper we describe a revised picture of gap closure and radiation physics in the outer 
magnetosphere, after reviewing briefly difficulties with the CHR scheme (§2). Although this 
radiation model shares the limitation of present cap models in that it does not start from a first 
principals computation of acceleration fields and particle densities, we find that a simple premise 
for the acceleration voltage allows a computation of the pulsar efficiency and its variation with 
spin parameters (§3) - the predicted luminosities compare well with those observed by EGRET 
and very high apparent efficiencies can be naturally accommodated in this model. In §4 we turn 
to the spectrum predicted by our radiation model. After describing the phase-average pulsar 
spectrum and dependence of spectral breaks on pulse parameters, we show how the high energy 
spectrum should vary through the pulse and present a test computation for Vela (§5). The 
computed phase- resolved spectra are compared with EGRET results. We conclude in §6 by noting 
3-D extensions to the spectral computations that should provide more precise models and by 
describing extensions of the pulse modeling to both lower and higher photon energies. 
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2. Outer Magnetosphere Radiation and Gap Closure 

CHR developed a version of an outer magnetosphere gamma ray emission zone in which 
acceleration in ‘Holloway (1973) gaps’ above the null charge surface brought a co- rotation charge 
density of particles to Lorentz factors 7 ~ 10 7 ' 5 . These primary charges (ei) were radiation-reaction 
limited by inverse Compton scattering on a posited dense bath of soft photons. The primary 
energy was thus radiated in a ~1-10 TeV Compton spectrum (71) that was fully absorbed via 
pair production on the same soft photon bath (73). The resulting secondary charges (e2) suffered 
strong synchrotron losses that resulted in a power law spectrum (72) of secondary radiation 
identified with the high energy pulses observed from the Crab and the Vela pulsars. In the CHR 
model the gap was assumed symmetric so that primary charges were produced all along the gap 
zone - one sign would accelerate and produce an outwards beam of radiation, while the opposite _ 
sign would similarly produce an inwards beam of 7-ray flux. In the final step of this picture, the 
opposing 72 beams provided for further pair production; the resulting low energy e 3 pairs were 
argued to produce the soft tertiary photon bath (73) needed for the original gap closure. 

Since this model was presented, several difficulties in explaining the observed pulsar radiations 
have arisen. For example, Chiang and Romani (1992) showed that the CHR picture in which a 
single bundle of field lines was active in the outer magnetosphere and which provided a double 
pulse by utilizing the opposing beams from opposite gaps did not produce pulse profiles like those 
observed. Moreover, in the original CHR model, it proved difficult to have a gap active for pulsar 
periods much larger than ~ 0.1s so that objects such as Geminga were difficult to understand 
(although Chen and Ruderman 1993 argued that highly inclined dipole rotators a ~ zr/2 might 
remain active in the Geminga regime). A key difficulty with this picture is that it was essential 
for the gap to be optically thick for pair production to absorbed the (unobserved) TeV photons 
produced by the the primary particles. This in turn requires very high luminosities for (73) 
which were in the near IR (Vela) to optical- UV bands (Geminga) - the required fluxes were 
generally orders of magnitude larger those observed from these objects (cf. Usov 1993). In fact in 
a computation of the CHR scheme, following the variation in the pair production and radiation 
efficiencies as particles moved through the outer magnetosphere gap, Chiang and Romani (1994) 
were unable to obtain sufficient tertiary flux for self-consistency, even for the Crab pulsar. Finally, 
the symmetry of the gaps suggests that the neutron star surface heating from the downwards 
directed particles and radiation might be excessive, as noted by Arons and Scharleman (1979). 
Usov (1993) showed that this argument applied to CHR- type outer magnetosphere acceleration 
produced significant excess thermal flux for objects such as Geminga. 

CR noted that if the outer gap spanned the entire set of field lines crossing the null charge 
• B surface above the surface of last closed field lines, then the outward directed beam from a 
single pole could produce pulse profiles quite similar to those observed from high energy pulsars. 

In RY, this geometrical model of pulsar beaming was further developed starting from the vacuum 
dipole solution and showing that observed pulse profiles, phase lags from the radio emission and 
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polarization properties could be naturally reproduced. Noting that the gap width controlled both 
the high energy beam extent and the 7-ray efficiency, YR showed how a synthesis of the galactic 
neutron star population produced gamma ray pulsars in good accord with those observed and 
implied that the bulk of the unidentified sources are indeed isolated r ~ 10 4 ' 5 y pulsars. For 
details of the beaming geometry see RY; the population work is extended to include more precise 
description of local high mass star formation, candidate identifications for the EGRET plane 
sources and comparison with model luminosity functions in Yadigaroglu and Romani (1996). 

To support these successes it is important to lay out the basic physics of gap closure 
for this outer magnetosphere model. Our model is based on curvature radiation in the outer 
magnetosphere and so has important antecedents in the work of Morini (1983) and Ochelkov 
and Usov (1980). There is also similarity to the gamma ray production in the slot gap of Arons 
(1981), although the basic geometrical and electrodynamic assumptions follow those of CHR more 
closely. Since the outer gap emission zones typically start from ^ 10 neutron star radii, high order 
magnetic multipoles and offset surface dipoles will have little effect on the magnetic field structure 
(although they might be quite important for the charge density boundary conditions at the lower 
edge of the gap). Similarly, emission generally is dominated by r L 0.3 r^c, where r^c = cP/ 2 ir. 
is the light cylinder radius. Accordingly, perturbations from particle inertia and currents will not 
be very strong in the low density regions of the gap. We can therefore, with some confidence, take 
the magnetic field geometry to be the retarded potential vacuum dipole solution (e.g. Deutsch 
1955) which lags increasingly from the static dipole solution as Q x r — * c. While tedious to write 
analytically, the computation of these field lines is straightforward on the computer and forms the 
foundation of our gap geometry. Observations may eventually demand perturbations from the 
vacuum approximation which will provide details on the magnetosphere current structure. 

We assume, following CIIR, that there is a starvation gap starting at the intersection of 
the null charge surface and the surface of ‘last closed field lines’ for which particle can co-rotate 
with the star. This closed surface is determined by selecting the field lines parallel to fl at a 
perpendicular radius r± £ r^c- The width w(4> az ) of this gap varies with azimuth around the 
polar cap formed by the projection of last closed lines to the star surface. The outer gap will span 
a fraction / wd<f) az of the field lines emanating from the cap and along these lines we assume that 
vacuum zone experiences a proportional fraction of the open field line potential drop, namely 

* gap ~\(^) 2 R .S,w (1) 

where B s is the polar cap surface field, and R* is the neutron star radius. The distribution of this 
potential drop along the gap depends on the details of the surrounding charge densities. CHR 
assumed a uniform drop, so that the acceleration field was Ey = $ gap / Rlc- If we assume that the 
field scales with vxB then we find Ey ~ r~ 2 . Because particles are radiation-reaction limited, 
and because the range of radii is modest, we find that results do not depend strongly on the form 
of Ey, but adopt Ey ~ r 0 ' 1 in the following treatment. 

We next assume that charges are created in this gap and that a modest fraction of the local 
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co-rotation (Goldreich-Julian) charge density 


tig-j = 7 x 10 10 B z (r)P 'em 3 


( 2 ) 


experiences this acceleration field. These charges are limited by curvature radiation to 

= ( 5.6XIO-O ^ " 51 * W 7 (X.J» 1 ,/P.,)*'V'* (3) 

where the surface cap field is B s = 2 x 10 12 i?i2G, the pulsar period is O.lfLiS and the local field 
line radius of curvature p is measured in units of ric • In general X e < 1 since both the overall gap 
width w and the presence of a screening charge density n £ uqj can decrease £y below its full 
vacuum value. Here we are interested in primary charges in the low density gap, so that X e = w. 
These charges emit a curvature spectrum with a cut-off energy 

E c oc 7 3 /p = 5.7 (X e B n ) 3 ' 4 P:l /4 GeV (4) 


for a typical radius of curvature p ~ r^cl 2 near ryve- These photons comprise the main observed 
EGRET spectrum of the 7-ray pulsars, with spectral cutoffs from ~ 1 - lOGeV. Note that in 
traversing the outer gap, the primaries each radiate 


Th*y 


eE\\rLC 
2 E c 


1.2 x 10 5 


x e B u y /4 

P - 1 ) 


( 5 ) 


peak photons. 


As in the original CHR picture we need pair production in the gap to tap the potential 
drop; with a magnetic field at r^c of B^c ~ 10 7 Pr^G only 7-7 pair production is viable. 
The multi-GeV curvature photons pair produce most efficiently on a ^ IkeV background. This 
is provided by the thermal surface emission. Figure 1 shows a range of theoretical thermal 
luminosities (measured at infinity) for models with non-exotic interiors (after Ogelman 1995) 
and r, = 10 6 cm. The corresponding surface temperatures remain above ~ 10 6 K for 10 5 y. For 
pulsar younger than ~ 10 4 y or with very high fields Ogelman notes that ‘magnetorotation ’ soft 
X-ray flux, powered by pulsar spindown and to be identified with the gap emission, exceeds the 
thermal surface flux shown here. The drop off in surface temperature after 10 6 y inhibits outer 
magnetosphere pair production. In our pulsar model primary photon aberration away from the 
field line direction ensures that pairs are initially produced with modest pitch angles ^ ^ 0.1 and 
7 ~ 10 3 ' 5 ; the perpendicular momentum is rapidly radiated away giving a synchrotron spectrum 
peaking at ~ 1 - lOMeV (§4) which dominates the soft gap emission. The charges are also 
accelerated via Ey to join the radiation reaction-limited primary population. As the particle 
population builds up beyond ncj , £y decreases - the resulting steep primary spectrum can give 
the curvature spectrum a softer index than the E'^dN/dE 1 ~ E* 4 ^ 3 expected for mono-energetic 
curvature radiation. 


The cross section for photon-photon pair production is 


ct 77 (E 1; £ 2 ) = ?r r 2 e {(5 cm {p 2 cm - 2) + (3 - /3 4 m )ln[7 cm (l + P cm )]}hln 


( 6 ) 



- 6 - 


where 7 cm = , j*Ei/mc 2 (l — (3*c os0i). With the total angle between the incoming photons as 
0 = 0 1 + 0 2 = 202 + 60 where tan(£0/2) = tan(#/2)(jEa — .E 2 )/(^i + E 2 ) we can solve for the CM 
frame velocity as (3* = ( E 2 — I?i)/(.E2COS02 - ^iCos^) giving <r^ r (Ei,E2,0). The result is a cross 
section of a 2 x 10 -25 cm 2 above a threshold E 2 of 

^ T m 

for primary 7 of energy E\ incident at angle 0. In practice, the optical depth for primary photons 
traveling through the magnetosphere is small £ 10 -3 — 10~ 2 . Thus while the bulk of the primary 
flux escapes to produce the observed beams, a small fraction of the n ^ ~ 10 5 primary photons 
produced by each primary charge produces pairs in the gap. The ‘gap closure’ is self regulating: 
if too many charges are produced, the vacuum width of the gap w decreases, E\\ and E c decrease 
and e t hr moves farther out on the Wien tail of the thermal surface distribution, decreasing the 
target photon cross section. Thus as particle are produced along a given field line the charge 
density builds up towards uqj ■ We define w(4> az ) to be the fraction of the open zone lines (at a 
given azimuth on the polar gap boundary defining a ‘last closed’ field line) with n < naj ■ Pair 
production of primary photons will continue outside of the vacuum region, but with decreasing 
efficacy due to the declining soft photon flux at r > ryvc- We thus expect the gap width to 
decrease with height over a distance £ r/vc to an asymptotic value for a given cap azimuth. The 
gap width will vary with azimuth around the polar cap. 

Computing the gap width in detail is difficult since the photon pair production depends 
critically on the angle 0. There are two effects important to its computation. First at high 
altitudes the primary photons will be aberrated ahead of their parent field line by the rotation and 
will thus have significant angle with respect to the largely radial surface flux - this increases the 
pair production efficiency of photons originating in the outer magnetosphere. Field line curvature 
will abet this somewhat and will cause leading/trailing asymmetries in the gap width. Second, at 
low altitudes the stellar radius provides a finite angle for the photon cross section. In some cases 
it is also significant that soft thermal photons will be resonantly scattered in the charge loaded 
regions of the closed zone. Thermal surface photons of energy E 2 will cross cyclotron resonance at 

r res = r/vsC^/ll-e^keV) 1 / 3 (8) 



in the dipole field around the star. Daugherty and Harding (1989) give the resonance part of the 
photon scattering cross section as 


c(E) 


aj 


E 2 


2 (E - E cyc y + T 2 e 


( 9 ) 


where E c = 11.61?i2keV and the natural linewidth due to excited state decay is 
Te = |am 2 (5/4.4 x 10 13 G) 2 = 2.6 x 10~ 3 5 2 2 (r)keV. This leads to very narrow resonant 
region and the scattering optical depth on the co-rotation density can be approximated as 
r ~ ^ 2 -(E%/T 2 E )(r res T E/3E2)naj(r res ) ~ 2 x 10 _2 [5i2/(£2/0.1keV)] 1 / 3 /P_i. Since the total e* 
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density near the neutron star is likely to be many times the net charge density uqj in the closed 
zone it is likely that the resonant scattering of the soft thermal flux increases the apparent size of 
the neutron star for the outgoing photons. These regions are Thompson thin unless n > I0 6 naj 
for these periods so there is an appreciable range of closed zone densities for which resonant 
scattering is the only substantial effect. This effect can be significant for low P, high field pulsars, 
so it is included although it does not generally affect gap closure for pulsars older than Vela. The 
effect of resonant scattering on the pulse profile of the thermal surface flux is described elsewhere 
(Romani and Rajagopal, in preparation). 

Computation of the evolution of the gap width with position in the magnetosphere involves 
following the population of curvature photons near E c as they propagate across field lines in the 
outer magnetosphere and tracing their conversion into pairs which subsequently travel along the 
field line at the pair production site. When the density from all pair production on a given field line 
reaches ~ hqj the gap will be closed. The increasing particle density near the gap upper surface 
defines the locus of field lines on which tangentially beamed curvature and synchrotron radiation 
define the observed pulse as described in RY. Note that the described gap is asymmetric ; since the 
bulk of the pairs in the high E|| region are produced near r ~ rpic, the outward-flowing charge will 
be re-accelerated and radiate n 7 photons in traversing the gap. The opposite, inward-accelerated 
charge will rapidly cross ryvc and leave the strongly charge depleted region. While we do then 
expect some inward-directed radiation, this flux will be much fainter than that of the outward 
beam and will be strongly beamed towards the spin equator. This should not contribute strongly 
to the observed pulsars viewed away from £ = tt/ 2. Moreover, except for very highly inclined 
dipoles, most of this smaller flux will not impact the stellar surface, leading to greatly reduced 
surface heating. If the synchrotron flux from these regions exceeds the thermal surface flux, as 
observations of the r ~ 10 3 y pulsars shows (Ogelman 1995), then computation of gap closure 
will be a difficult self-consistency problem. However, when thermal surface flux dominates in the 
critical lOOev-lkeV pair production region then simple approximate sums for the gap closure and 
pulsar efficiency can be obtained. 


3. Gap Closure and Efficiency Evolution 


An outer magnetosphere model should depend only on the spin parameters P and P 
(represented here for comparison with data by characteristic age r and dipole magnetic field 
B ) and the magnetic inclination angle a. Perturbations such as surface multipole structure, 
spin-history, or external photon fluxes should have only weak effects in the gap zone. Although we 
can in principle compute the entire 3-D gap structure and pulsar beams, we start by computing 
characteristic gap sizes and beaming angles that will provide a pulsar 7-ray efficiency 7 7 (r, B,a) 
according to the physical scheme in §2. 
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The characteristic scale in the outer magnetosphere will be 

ric = 4.8 X 10 8 PTi cm (10) 

The intersection of the null charge zone with the last closed field line will vary with azimuth but 
we adopt here the estimate for the radius of the null charge point in the Cl - B plane given by 
Halpern and Ruderman (1993), r^c/ r LC ~ (3 tana/2)~ 2 for a £ 7r/4. For a = 0, r^c = 2/3. As 
the gap grows the null charge intersection of the gap upper surface moves outward. Noting that 
we can relate the ‘closure radius’ of the gap lines to the polar cap angle as in the static dipole 
solution and that the footpoint angle of the gap upper surface scales sis (1 — w) we write the radius 
of the gap closure point as 

r, « 2.1 x 10 8 Pr 1 1 [(l - w)tana] -2 cm. (11) 


The finite size of the scattering surface in the closed zone ensures that the minimum photon-photon 
angle is ~ r res /2 1 / 2 r, where r rea is given by (8). At higher altitudes the photon angles are 
dominated by aberration with 6 ss 7rr,/2r£,c ss 0.7[(1 - w)tana] -2 . 


Taking the typical field line radius of curvature in the gap to be p ~ r^c/2, we can 
adopt (4) and (5) for E c and n 7 and compute the pair production on the thermal surface 
flux. The evolution of the neutron star surface temperature is adequately described by 
T(t ) = 10 6 /v (t/10 5 y)“ OO5 exp( — t/10 6 y) (giving the solid line in Figure 1) out to several times 
10 6 y. Complications due to non-thermal surface spectra induced by neutron star atmospheres (eg. 
Rajagopal and Romani 1996, Shibanov et al. 1992) are ignored here as is additional soft X-ray 
flux produced by the gaps themselves. Comparing the pair production threshold energy (7) with 
the thermal peak, and taking cos 6 « 9 2 / 2, we see that 


tthr 

w 


> 7/4 

-1 


2A 0 2 T 6 (X E l B n )W (12) 

Thus for typical pulsar conditions we will be on the Wien tail of the thermal surface spectrum and 
the soft photon density above threshold will be 

2aT 3 .e thT .o_^ t 


/ '-inr \ 9 / \2 

n 2 * e kT ( 7 '«/ r *) • 


(13) 


Finally to compute the gap width we compute the pair multiplication factor n mu i = n 1 n 2 (J T1 Ti\ 
taking X E = w, we can then solve n mu i(w) = 1 for the width at which particle multiplication is in 
equilibrium and gap closure occurs. 


In this model, the total potential drop (and hence the energy deposited in curvature photons 
per primary) scales as w. Moreover the total number of charges participating will be proportional 
to the gap width so we expect the efficiency to be quadratic in w. Since w is a linear measure 
of width and since the field lines initially curved away from the spin axis are the ones closed at 
moderate altitude, we adopt' 

= w{w - y)n mu /(u;) 


(14) 
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which accounts for the convergence of the polar cap lines to the magnetic axis and where the final 
factor of n mu i accounts for the low charge density in only partially closed gaps. The model is 
clearly not complete for very large w or very small a, so we restrict r,- to be < r^c and w < 0.8 in 
numerical computations. 

In addition to rjy, the fraction of the spindown power converted to MeV-GeV photons, it is 
important to estimate / 7 , the beaming factor of these 7-rays on the sky. As shown in YR, this 
determines both the solid angle over which a given pulsar will be observable (centered on the 
spin equator) and the apparent 7— ray luminosity for observers within the beam. For example, if 
the true beaming factor is / 7 but a beaming factor of 1 is assumed (i.e. isotropic radiation), the 
inferred pulsar efficiency will be 77 = 77-,/ / 7 . Although a detailed model for the 7-ray beam at each 
r, B and a is needed for a precise value for 77 7 we can adopt here the following fit to the models 
plotted in YR, Figure 1 

_ / a \ 1/2 (1 - w) 

h ~ Uo°; (1 + w) 

which is accurate to better than ~20% over most of the parameter range. 

Together these estimates give rj y (T,B,a ) (and L 0 b s = EsDV-y/i^nd^fy)). Figure 2 shows the 
efficiency variation with characteristic age for a = 65° and several values of B. In addition to the 
true efficiencies, we show the efficiencies inferred for a beaming factor of 1 steradian, along with 
the E > lOOMeV efficiencies inferred for several EGRET pulsars (Fierro, 1996), again with an 
assumed / 7 = l/4ir. Several features are worth noting. First, efficiencies increase strongly with 
spindown age and are also higher for large B objects at a given r. This is because for large B and 
or r the pulsar spin period P = 0.248But^ 2 s increases, so that the light cylinder radius and t-jvc 
increase. This makes gap closure more difficult; with a larger w a higher efficiency and smaller / 7 
result. The steep fall-off in r] y for r < 10 4 y is due to the relatively large soft photon angles imposed 
by the circumstellar resonant scattering. If the closed zone does not contain sufficient density to 
allow this scattering, 77 7 will be somewhat larger in this region. Note however that for high B this 
is also the region where magnetospheric flux can dominate surface emission; this unincluded flux 
will conversely decrease somewhat the predicted efficiencies for very young pulsars. The steep 
fall-off at large ages is due both to the growth of r, and the rapid fall-off in the surface T beyond 
10 6 y. By comparing the true efficiency curves with the curves for / 7 = l/4x, one sees that the true 
beaming factors are large for inefficient pulsars, but drop to ^ 10% for highly efficient pulsars. 
Thus, although the true 7-ray efficiencies are only 0.2 — 0.3 for these objects, if isotropic emission 
is assumed one will infer 7-ray luminosities exceeding the total spin-down power. 

In figure 3 we display the dependence on magnetic inclination a. For old pulsars large 
inclinations provide for small r t . This in turn limits the contribution of aberration to 9 and 
promotes high 7-ray efficiency. In computing n mu i it is found that in this regime the gap self-limits 
so that the exponent in (13) is nearly constant. This implies that P 7 / 4 oc B 3 / 4 . Substituting in 
for the evolution of P this in turn assures that r 7 / 8 oc B ~ 1 . Therefore pulsars with tB ~ constant 
have similar efficiency variation. For younger pulsars the algebraic factors become significant and 
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the behavior is more complex; peaks near a ~ 65°. Note that, as in Figure 2, the rj 0 k 3 inferred 
for inefficient pulsars are much smaller than their true values if / = \/Air is assumed. 

4. Radiation Spectrum 

The curvature radiation power spectrum from a primary with energy 7 e mc 2 is 

P(E) = E 2 ^- = * 1/2 -p eE F syn (E/E c ) = 6.02 x 10 7 ^^F syn (E/E c )MeV/s (17) 
atat, tip p 

where E c is the cuttoff energy (4) and F syn (z) = x / y °° I( 5 / 3 (y)dy. For E << E c this gives the 
steeply rising spectrum P(E ) oc E 4 ! 3 . While recent observations of 7 -ray radiation from efficient 
old pulsars and phase- resolved analysis of other objects indicates photon number spectra as steep 
as dN/dE ~ E ~ 13 , it has long been known that younger objects such as Crab and Vela exhibit 
a softer spectrum. Of course, a steep spectrum of primary particle energies can produce a flatter 
curvature spectrum over a modest range. 

As charges are produced and accelerated above r, we expect two effect that can generate a 
spectrum of radiation reaction limited particles. First, near r{ particles on a given field line will 
experience a range of E\\ as they pass to regions of decreased w and higher charge density. Second, 
the J5|| should not drop abruptly to zero, but should exhibit some penetration of the increasingly 
dense plasma at the gap upper surface. For concreteness, we posit that X e oc ncj/n for n > ncj 
and that charge increase has ~ constant scale length transverse to B. A careful treatment of the 
charge and field distributions around the gap will be needed to see if the anomalous resistivities 
required for this field penetration ansatz are present, or if we will need to rely on the acceleration 
history at the gap onset to supply the electron energy spectrum. Note that our ansatz gives 
dN e /dE\\ oc .Ejj " 1 and thus as long as the particles are radiation-reaction limited (3) this gives rise 
to a steep electron energy spectrum dN e /dE e oc E~ p with p = 4. The electron spectrum will 
reach radiation reaction energies as long as t acc = -pmc 2 /eE\\c is less than t// = P/2ti. For the 
acceleration potential ( 1 ) this leads to the condition 

p 2 

X E > 3.8 x 1(T 3 — (18) 

77-012 

to reach Lorentz factors of 10 7 77 . This electron spectrum causes the curvature spectrum 
to develop a broad maximum with power spectrum E 2 dN/dE ~ £-( p - 5 )/ 3 ^ E 1 / 3 ( i.e . 
photon number spectrum dN/dE ~ E~ 1 ' 7 ) This spectrum will extend from E c ~ 3GeV to 
Eb ~ EcX^I^^ ~ 20MeV for parameters similar to Vela’s. Below this break energy the curvature 
spectrum will steepen towards the P(E) ~ E 4 ! 3 from single particle curvature radiation. For 
longer period or lower field pulsars such as Geminga and PSR1055-52 XE,min w iU be much 
larger and the GeV spectrum will approach E~ 1 as the particle spectrum becomes more nearly 
mono-energetic. In general, we expect the spectral index to be correlated with Tj 1 . 
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Since photons are aberrated across field lines, pairs are initially produced with significant 
pitch angles ip ~ 0.1 and initial Lorentz factors 7 ,• ~ E c / 2 ~ 3 X ■ Near the 

null charge line these will synchrotron radiate, delivering a spectrum with an upper cutoff 

Esyn = 0.5^ /2 5 1 5 2 /2 Pr i 13/2 (r,/0.1r z , c )- 2 MeV (19) 

taking as a rough approximation ip ~ Ti/r^c. For Vela parameters, this gives E ayn ~ 2MeV, 
while for the Crab E ayn ^ 300MeV. Since E\\ acceleration does not counter pitch angle losses, the 
electrons cool, resulting in dN/dE± ~ E± p with p ss 1. The corresponding synchrotron spectrum 
is is P ayn (E ) ~ £-( p - 3 )/2 ^ E +1 . This flatter spectrum will dominate the curvature emission for 
E < E ayn . Note that since P ayn {E ) oc B(E/E ayn ) 1 ^ 3 ~ B 1 ^ 6 there will be a weak increase in P ayn 
at a given E for high field pulsars. We can also estimate the ratio between the total curvature 
power P cr ~ Eyn-y ~ 4> ffap and the total synchrotron power P ayn ~ 7 ,n mu / ~ £ ; 7 n t / l (r,)cr- y7 r,- by 
noting that P ayn /P C r ~ t 77 ~ 3 x 10 -3 for gaps controlled by thermal surface flux. For young, high 
field pulsars such as the Crab and PSR1509-58, E ayn will move above lOMeV and synchrotron flux - 
will dominate the thermal surface emission. In this case, we expect that n^r,) will substantially 
exceed that due to the thermal surface flux and r 77 will increase towards unity, giving a power law 
non-thermal spectrum at keV-MeV energies, and spectra peaking below the EGRET range. Figure 
4 show the schematic phase averaged pulsar spectrum, including thermal surface flux, curvature 
radiation and synchrotron emission for r ~ 10 4 5 y and B \2 — 3. 


5. Phase- resolved EGRET spectra: Vela Model 


The pulsar spectrum of the previous section shows many features present in the phase 
averaged spectra collected by Thompson (1996). It is, however, important to remember that these 
estimates are based on characteristic values for r,-, 9, etc. Since these quantities vary with position 
in the magnetosphere we expect spectral variations through the pulsar profile. Such variations 
can be noted by comparing pulse profiles from different energy ranges - for example in our model 
we expect power-law synchrotron emission from near to be present in the X-ray light curves of 
young pulsars leading to power-law tails in ROSAT spectra, peaked between the 7-ray maxima for 
viewing angles ( not too far from 7r/2. We also can follow variations in the curvature spectrum 
near the peak. For the brighter EGRET pulsars, high quality phase-resolved spectra are now 
available so we can compare the predictions of our outer magnetosphere radiation model with 
these data. 

To implement the gap closure model of section 2 in a local manner, we must follow the photon 
beams as they propagate to the light cylinder, crossing field lines and converting on the thermal 
photon flux. The build-up of charges on a given field line then determines the locus of the gap 
radiation surface and the local value for w. Such models produce asymmetric gaps with variation 
along the surface of last closed lines (in preparation). As an intermediate step we can however 
trace photons and charges along a single line to confirm the basic behavior described in §2 & 3. 
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As an example we will estimate gap closure and radiation for the Vela pulsar, which has the best 
count statistics in the EGRET data. 

The gap closure sums of section 3 predict w = 0.17 for Vela-like parameters (Log(r) = 4.04, 
i?i2 = 3.3, a « 70°). Adopting the the potential drop 4> aap , distributed as r -1 we have computed 
an approximate value for w along several lines spaced around the polar cap by taking the curvature 
photon flux to be that produced at lower altitudes on the same line. These photons are allowed 
to pair produce on the thermal surface spectrum, using the detailed kinematics of (6). For field 
lines curved towards the equator, we find that the gap reaches a stable width at w = 0.17 — 0.24; 
this asymptotic value is reached within a pathlength 6s < along the asymptotic field 

line. This supports the values for the gap closure estimates in §3. Since cross field propagation 
is not presently included, we cannot examine the variation of the gap width across the outer 
magnetosphere radiation zone. We can, however, estimate the phase averaged spectrum and 
illustrate variations in pulse phase by following the changes in the emission location. 

We define the last closed field lines as those parallel to Cl at r^c/ 2 1 / 2 (where the co-rotation 
velocity equals the Alfven speed, cf. Henrickson and Rayburn 1974), this provides a polar cap 
~ 2 1 / 4 times larger than the simple vacuum value. Then adopting a constant w — 0.2 gap width, 
a = 65° and a rotation axis- line of sight angle C = 79° as in RY we find the location of the 
emission region on the upper surface of the gap viewed as a function of pulsar phase. The altitude 
of the emission point at each phase along with the distance from the null charge surface 6s are 
shown in Figure 5. Superposed is a light curve generated assuming uniform emissivity along the 
gap surface for reference. The radio pulse occurs at phase 0. Note that in the interpulse region 
emission occurs lower in the magnetosphere and closer to the null charge surface. 

Our trace of w(r ) along a typical field line allows us to estimate E\\(r) in the gap, but does not 
accurately reproduce the build-up of charge beyond n ~ ngj. We therefore allow the maximum 
charge density to increase oc Exp(5s/0.1r£,c) 2 as pair production continues over a distance ~ r,. 
Charges are injected with 7 = 10 7 ' 5 and allowed to come into radiation reaction equilibrium as 
they flow along the gap. This prescription allows the particle population to extend to lower 
energies as XE,min oc 1 /N decreases higher in the gap. The result is a narrow distribution of 
electron energies near r,-, spreading over a factor ~ 3 in E e as 6s approaches 0.3 r^c- We next 
compute the full curvature spectrum as a function of pulse phase for the Earth line of sight at 
C = 79° and compare with observed EGRET spectra. 

In figure 6, EGRET phase resolved spectra from a pulsed analysis (Fierro 1996) are shown for 
the first peak, the second peak and the interpulse region. Pulsed analysis avoids contamination 
from a soft spectral component, largely unpulsed, which is likely background flux and also allows 
spectra to be extended above lOGeV. The model spectra, computed for the same phase range 
as the corresponding data, reproduce the details of the phase variations quite well, although the 
components have been individually normalized. In the constant w model described, the first pulse 
is slightly too bright and the second pulse is too weak by a factor of 3. It is likely that this will 
be rectified in a 3-D computation with variable gap width, since the leading emission of PI from 
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large r is aberrated more strongly, leading to a narrower gap, while the increased width in P2 
region will result in a larger local curvature efficiency. It is encouraging that this test computation 
produces spectral slope variations in the lOOMeV to GeV range very close to those seen and phase 
varying high energy roll-overs consistent with the data. Note that this model predicts significant 
10-30GeV flux in contrast to polar cap pair cascade models. 

The phase averaged sum of the pulse components is shown in figure 6 and over a broader 
energy range in figure 7, combining data from three CGRO instruments. Here we also plot the 
estimated synchrotron spectrum with E syn following (19) and the synchrotron/curvature flux 
normalization F 3yn IF CT = 10 -3 ' 5 . Changes in the viewing angle can significantly affect this soft 
component, which here rises above the E < Eb portion of the curvature spectrum. 


6. Discussion and Conclusions 


We have outlined the physics of gap closure and curvature and synchrotron radiation in the 
outer magnetosphere. This model predicts efficiencies and beaming fractions for GeV pulsars in 
good accord with the data. Significant dependence of on the pulsar magnetic field was found. 
Since the 7 -ray emission is also quite sensitive to the magnetic inclination a (which, along with 
the spin parameters B and r, completely determines the pulsar properties), careful polarization 
observations and models of the radio profile to constrain a for young pulsars are particularly 
important. 

A test computation with a constant gap width w in this pictures produces an encouraging 
match to the spectra and phase variations of the Vela pulsar. These results should be substantiated 
by computations following the variation of w for various gap field lines and for other pulsars (in 
preparation), but it seems likely that the modeled spectral variations will be quite robust. The 
tendency of older, more efficient pulsars towards harder spectra is natural in this model; photon 
number spectral indices approaching -1 can be expected for some phases of old pulsars, but there 
should always be significant curvature in the lOOMeV-GeV range so power laws will not provide 
adequate fits to high quality data. 

The computations in this paper have been directed to pulsars with ages greater than or 
similar to that of Vela for which thermal surface emission should dominate the flux in the sub-keV 
range. Pulsars such as Crab, PSR0540-69 and PSR1509-58 are clearly dominated by non-thermal 
power laws at these energies. We infer that this flux is synchrotron radiation from the outer 
magnetosphere and that larger r 7 - K will assure copious pair production and small w and GeV 
7 — ray efficiencies. Computations of gap closure will be more difficult in this case since, as in the 
original CHR model, both the high energy photons and the soft target flux are produced by the 
gaps themselves. Closure itself will, however, remain much easier than in the CHR scheme because 
of the modest pair creation optical depth required. 

Since these r ~ 10 3 y pulsars will have pair creation extending well above r^c it is tempting 
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to interpret likely spectral differences in the light of observations of the Crab pulsar. First, with 
T~n & 1 we can expect the magnetospheric pair creation rate to be n^wA cap nojc ~ 10 38 e ± /s for 
Crab parameters. This is consistent with the pair content inferred from observations of the plerion, 
although further particle acceleration is required in the wind shock. The continued production 
of pairs out to r ~ r^c gives rise to the sharp optical pulses from the outer gap region. The 
Crab gaps are thus also much more symmetric than those of Vela and other older objects - weak 
absorbed emission from inward directed beams crossing the magnetosphere may be expected. 

It also seems likely that the dense pair plasmas in the high altitude gaps are what give rise to 
the coherent radio emission from these regions exhibiting giant pulses and spectral peculiarities. 
Finally, there is the question of the field lines that curve away from the pulsar spin equator, cross 
the spin axis and reach the Q ■ B surface near r^c opposite the standard outer gap zones (these 
connect with the normal gap at large radii, but may be thought of as equivalent to the ‘short’ 
gaps of CHR, figure 1). Normally these gaps do not close well and contribute only weakly to 
the 7 ray emission. However, if we assume that the strong pair creation at large r, in the Crab 
allows these to be active and trace the expected pulsed flux for the Earth line of sight, we find 
two pulse components lagging P2 by ~ 60° and ~ 100°. It is intriguing that recent 4.8GHz radio 
observations show two additional pulse components at these phases (Moffet and Hankins, private 
communication). It will be interesting to see if these component exhibit giant pulses and related 
radio phenomena and to see if careful analysis of apparent off-pulse structure in the EGRET Crab 
data can detect weak 7-ray emission at these phases. 


Finally, we can estimate the inverse Compton scattering of the primary electrons on soft 
photons from the pulsar gap. The 7 ~ 10 7,5 primaries will most efficiently upscatter 10~ 2 eV (IR) 
photons from the low energy end of the synchrotron spectrum. Most of this scattering is in the 
Klein-Nishina regime where Blumenthal and Gould (1970) give the upscattered number spectrum 
as 


dN/dtdE = 2Tr ^’‘ MJdE, 
' 7e £2 


2? In q -(- (1 + 2q)(l - q) + 


(r e <?) 2 (i-<?) 

2(i + r e9 ) 


( 20 ) 


with r e = 4jE27e/ J7 ic 2 , q = p an d £1 = E/‘y e mc 2 . Without a careful computation 

of the synchrotron beaming, it is difficult to estimate the photon density n{Ei) at these 
energies in the gap. Roughly, though, we can expect E y dN~,/dE^ ~ const for our synchrotron 
spectrum. Comparing with the flat power law synchrotron flux at thermal peak, we typically have 
nth, ~ 10np£, for 10 4 y, £12 = 3 pulsars while crj ~ 30 ^, so we expect tjcs ~ 3 x 10 -3 , scaling from 
the estimates of section 4. Integrating the upscattering over the primary population N(-y e ) down 
to the cut off from the minimum Xg (Eq. 18), we obtain the TeV Compton spectrum shown in 
Figure 4. This pulsed TeV flux containing somewhat less than 1% of the pulsed GeV flux and is 
an attractive target for future air Cerenkov experiments. 
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Figure Captions 

Figure 1 . Soft X-ray luminosity evolution for young pulsars. The solid line approximates 
theoretical neutron star cooling curves for stars with non-exotic interiors (model ranges - points). 
The dashed curve shows the empirical ’magnetorotation’ luminosity law proposed by Ogelman 
(1995) for B \2 = 5. 

Figure 2. Efficiency (t/ 7 ) evolution for the high energy flux of young pulsars with various 
B and with a = 65°. Also shown are the inferred luminosity variations r)\ (dashed lines) for an 
assumed fixed / 7 = 1/4 x. Data points give the corresponding inferred 100MeV-2GeV efficiency 
ranges from EGRET observations of 7-ray pulsars (closed circles) and candidate unidentified 
sources (open circles) from Fierro (1996); again / 7 = l/4ir is assumed. For comparison the dotted 
line gives the single parameter efficiency estimate of YR: t? 7 (t) = 3.2 x 10~ 5 r 0 76 . 

Figure 3. True (solid lines) and 1 steradian beam (/ 7 = l/47r) estimates (dashed lines) for 
7-ray pulsar efficiency variation with magnetic inclination angle a. Curves show the variation for 
several combinations of B and r. 

Figure 4. Characteristic phase average spectrum for a young 7-ray pulsar (B i2 = 3, 

E c = 3GeV). Solid lines show the curvature spectrum (CR), the synchrotron emission near r, (Sy) 
and the thermal surface flux (kT) described in section 4. For comparison, the dotted curve gives 
the mono-energetic primary curvature spectrum with the same E c as in (CR). The dashed curve 
gives the TeV pulsed spectrum from Compton upscattering of the synchrotron spectrum on the 
primary e ± . 

Figure 5. Variation of the emission point as a function of pulse phase for Vela. Crosses give 
the radius of the (lower) emission point in units of r^c', open circles give the distance from the 
null charge surface in the same units. For reference, the solid line gives a pulse profile for uniform 
emissivity along the gap. 

Figure 6. Phase resolved spectra for the Vela pulsar. Data points are from a pulsed analysis 
for the first peak (PI), an interpulse zone (IP2), the second pulse peak (P2) and the total flux 
(Fierro 1996). Model curves for the same phase intervals and Vela parameters have assumed a 
constant gap width w = 0.2 and followed the evolution of the electron and photon spectra along 
the resulting gap. 

Figure 7. Phase averaged spectra for Vela compared with CGRO data. The solid lines give 
the curvature and synchrotron spectra (sum as dashed line) for the Vela model observed at the 
earth line of sight. 
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